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Abstract. 

We develop and apply a fully covariant 1 + 3 electromagnetic analogy for gravity. 
The free gravitational field is covariantly characterized by the Weyl gravito-electric 
and gravito-magnetic spatial tensor fields, whose dynamical equations are the Bianchi 
identities. Using a covariant generalization of spatial vector algebra and calculus to 
spatial tensor fields, we exhibit the covariant analogy between the tensor Bianchi 
equations and the vector Maxwell equations. We identify gravitational source terms, 
couplings and potentials with and without electromagnetic analogues. The nonlinear 
vacuum Bianchi equations are shown to be invariant under covariant spatial duality 
rotation of the gravito-electric and gravito-magnetic tensor fields. We construct the 
super-energy density and super-Poynting vector of the gravitational field as natural 
C/(l) group invariants, and derive their super-energy conservation equation. A 
covariant approach to gravito-electric/magnetic monopoles is also presented. 



PACS numbers: 04.30.-w, 04.20.-q, 11.30.-j, ll.30.Fs, ll.15.Tk 



Introduction 



There is a surprisingly rich and detailed correspondence between electromagnetism and 
General Relativity, uncovered in a series of fundamental papers by Bel |]l|, Penrose [Q] 
and others p, ^, ||, H, 0, ||, ^ (see ||^, |T0[ for more references), and further developed 
recently (see, e.g., [0, [12], |13|, |I4|, |15|, |T6|, |T^, [T8|, [1^). This correspondence is reflected 
in the Maxwell-like form of the gravitational field tensor (the Weyl tensor), the super- 
energy-momentum tensor (the Bel-Robinson tensor) and the dynamical equations (the 
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Bianchi identities). Another form of the correspondence arises in the search for geons 
(locahzed, non-singular, topological solutions of Einstein's field equations with mass 
and angular momentum): in the known (approximate) solutions, the geometry of the 
electromagnetic geon is identical to that of the gravitational geon [^, pH| . 

Here we pursue the 'electromagnetic' properties of gravity in areas which have 
already proved useful for extensions of electromagnetism to non-Abelian gauge theories 
and string theory. Our emphasis is on a 1 + 3 covariant, physically transparent, 
and non-perturbative approach, with the gravito-electric and gravito-magnetic spatial 
tensor fields as the fundamental physical variables. Using an improved covariant 
formalism, including a covariant generalization to spatial tensors of spatial vector 
algebra and calculus, we show in detailed and transparent form the correspondence 
between the electric/ magnetic parts of the gravitational field and of the Maxwell 
field. We identify gravitational source terms, couplings and potentials with and without 
electromagnetic analogues, thus providing further physical insight into the role of the 
kinematic quantities shear, vorticity and four-acceleration. 

In the vacuum case, we show that the nonlinear (non-perturbative) Bianchi 
equations for the gravito-electric and gravito-magnetic fields are invariant under 
covariant spatial duality rotations, in exact analogy with the source-free Maxwell 
equations for the electric and magnetic fields. The analogy is of course limited 
by the fact that the Maxwell field propagates on a given spacetime, whereas the 
gravitational field itself generates the spacetime. The electromagnetic vectors fully 
characterize a Maxwell solution, and duality maps Maxwell solutions into Maxwell 
solutions. The gravito-electric/magnetic tensors are not sufficient to characterize 
covariantly a solution of Einstein's equations - one needs also the kinematic quantities 
which are subject to the Ricci identities [^. Duality is an invariance only of the Bianchi 
identities, and not the Ricci identities, so that it does not map Einstein solutions into 
Einstein solutions. Nevertheless, the covariant gravito-electric/magnetic duality reveals 
important properties of the gravitational field. 

The covariant 1 + 3 duality has not to our knowledge been given before. Although 



duality invariance follows implicitly from Penrose's spinor formalism [g, this is in 



terms of the 4-dimensional Weyl spinor, rather than its 1 + 3 electric and magnetic tensor 
parts. Four-dimensional covariant tensor approaches to the electromagnetic analogy 
have been developed (see e.g. fl^), and non-covariant linearized Maxwell- type equations 



are well established, both in terms of gravito-electromagnetic vectors (see e.g. |23, 24 



and tensors (see e.g. [§). In [^, a covariant and nonlinear vector approach is developed 
for stationary spacetimes. Our approach is fully covariant and non-perturbative, and 
in addition is centred on the gravito-electromagnetic spatial tensor fields, allowing for a 
more direct and transparent interpretation based on the Maxwell vector analogy. This 
approach is a development of the work by Triimper p[. Hawking M and Ellis [0], and is 
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related to recent work on a covariant approach to gravitational waves [jl6|, |18|, |T9| and 
to local freedom in the gravitational field [|19]. A shadow of our general duality result 
arises in linearized gravitational wave theory, where for vacuum or de Sitter spacetime, 
there is an interchange symmetry between the gravito-electric and -magnetic tensors 

Duality invariance has important implications in field theory in general. It was 
essentially this symmetry of the Abelian theory, and attempts to extend it to include 
matter, which led to the Montonen-Olive electromagnetic duality conjecture that there 
exists a group transformation mapping electric monopoles into magnetic monopoles 
within the framework of a non- Abelian (specifically SU{2)) gauge theory |2^. This 
conjecture has proved particularly fruitful, stimulating work on S, T and U dualities 
in string theory (see e.g. ||2^, ^ HH), the extension of the electromagnetic duality to 
magnetically charged black holes and nonlinear electrodynamics []3ll , p2| and leading to 
the Seiberg-Witten proof of quark confinement in supersymmetric Yang-Mills theory 
via monopole condensation . 



We use the covariant spatial duahty to find the gravitational super-energy density 
and super- Poynting vector as natural group invariants, and derive a new covariant super- 
energy conservation equation. Finally, we discuss gravito-electric/magnetic monopoles, 
providing a covariant characterization, in contrast to previous non-covariant treatments 
34|, p4[. In the linearized case, we show that the Taub-NUT gravito- magnetic 



monopole given in I^S] is related to the Schwarzschild gravito-electric monopole by 
a spatial duality rotation and an interchange of four-acceleration and vorticity. This 
provides a covariant form of the relation previously given in non-covariant approaches 
It is well-known that the NUT metrics may be obtained from the Schwarzschild 



metric via the Ehlers-Geroch transformation 351. This transformation is in fact the 



generator of T-duality in string theory, but it is not a duality transformation in the 
sense described here, since it maps Einstein solutions to Einstein solutions and thus 
necessarily involves kinematic and geometric conditions in addition to a duality rotation. 
Furthermore, the Ehlers-Geroch transformation requires the existence of a Killing vector 
field, whereas the general duality that we present does not require any spacetime 
symmetry. 



2. Covariant spatial vector and tensor calculus 

To elaborate the electromagnetic properties of the free gravitational field in General 



a 



Relativity, we first present the required covariant formalism, which is based on p7 
streamlined and extended version of the Ehlers-Ellis 1 + 3 formalism . Then we give 
the covariant form of the Maxwell spatial duality in a general curved spacetime. In the 
following section we extend the treatment to the gravitational field. 
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Given a congruence of observers with four-velocity field m", then hat = gab + UaUb 
projects into the local rest spaces, where gab is the spacetime metric. f The spatially 
projected part of a vector is 

V{a) = ha^'Vb, 

and the spatially projected, symmetric and tracefree part of a rank-2 tensor is 
The spatial alternating tensor is 

^abc Vabcd^ ^[afec] ; 

where rjatcd = Vlabcdi is the spacetime alternating tensor. Any spatial rank-2 tensor has 
the covariant irreducible decomposition: 

^ab = ^hcdA'^'^hab + A(^ab) + ^abc-^^ , 

where 

^a 2^a,bc-^ ^ 

is the vector that is the spatial dual to the skew part. Thus the skew part of a spatial 
tensor is vectorial, and the irreducibly tensor part is symmetric. In the 1 + 3 covariant 
approach 0, all physical and geometric variables split into scalars, spatial vectors 
or spatial tensors that satisfy Aab = ^(afe)- From now on, all rank-2 spatial tensors will 
be assumed to satisfy this condition. 

The covariant spatial vector product is 

[V,W]a = eabcV'W\ 

and the covariant generalization to spatial tensors is 

B]a = Sabc^^dB'"^ , 

which is the vector that is spatially dual to the covariant tensor commutator. 
The covariant time derivative is 

A'^-b... = u'V.A'^-b... , 

and the covariant spatial derivative is 

Da^^ c--- = ha^h^q ' ' ' h/ ■ ■ ■ VpA"^ . 

t We follow the notation and conventions of 0. (Square) round brackets enclosing indices denote 
(anti-) symmetrization, while angled brackets denote the spatially projected, symmetric and tracefree 
part; a, 6, • • • are spacetime indices. 
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Then the covariant spatial divergence and curl of vectors and rank-2 tensors are defined 



by 13, [19 



div V = , curl K = Sabciy'V , (l) 

(div A)a = B'Aab , curl Aab = ecd(aD^A,)^ , (2) 

where curl Aab is tracefree if Aab = A{ab) ■ The tensor curl and divergence are related by 

The kinematics of the w'^-congruence are described by the expansion O = D°'Ua, the 
shear aab = D(aM6), the vorticity Ua = — |curlMa, and the four-acceleration Ua = U(a)- 
The above operators obey the covariant identities 

(D„/)- = Daf - |eD J + iiaf - (Ta'BJ - [uj , D/]„ + M^M^D J , (3) 

curlD,/ = -2M, (4) 

B^IV, W]a = Wcml Va - 1/"curl Wa , (5) 

B''[A,B]a = B^'cmlAab - A'^'cmlBab, (6) 

together with far more complicated identities [S^, |1^. In the case where spacetime 
is almost spatially isotropic and homogeneous, i.e. a linearized perturbation of a 
Friedmann-Lemaitre-Robertson- Walker (FLRW) background, some of the main further 



identities take the linearized form |31 

(D"K)' ^B-Va-HD^Va, (7) 

(D'Aab)' ^D'Aab-HD'Aab, (8) 

(curl Va)' ^ curl K - ifcurl K , (9) 

(curM^b)' ^ curli^b - if curl , (10) 

D^curlK ~0, (11) 

D'cml Aab ^ icurl (d^A,,) , (12) 

curl curl K ~ - D V, + (dV,) + | (p - 3H^) K , (13) 

curl curl Aab ~ - ^^Aab + |D(„D^A), + (p - SH"^) Aab , (14) 

where H is the background Hubble rate, p is the background energy density and 
-q2 _ j^aj^^ jg covariant Laplacian. 

The electric and magnetic fields measured by u"" observers are defined via the 
Maxwell tensor Fab by 

Ea = FabU^ = E(^a) , Ha = \eabcF^'^ = *FabU^ = Hl^a) , (15) 

where * denotes the dual. These spatial physically measurable vectors are equivalent to 
the spacetime Maxwell tensor, since 

Fab = 2U[aEb] + EabcH" . (16) 
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Maxwell's equations V[aFhc] = and V^Fab = Ja are given in 1 + 3 covariant form 
for Ea and Ha by Ellis [|^. In the streamlined formalism, these equations take the 
simplified form 

D"Ea = -2uj''Ha + Q, (17) 

D'^Ha =2u''Ea, (18) 

E^a) - curl Ha= - iQEa + OabE'' - , E]a + [m, H]a - Ja , (19) 

H^a) + curlE, = - iQHa + - [u, H]a - [il, E]a , (20) 

where g = —JaU"" is the electric charge density and ja = J {a) is the electric current. 
In flat spacetime, relative to an inertial congruence {Q = iia = uja = (Tab = 0), these 
equations take their familiar non-covariant form. 

Introducing the complex electromagnetic spatial vector field 2a = Ea + iHa, we see 
that in the source-free case (J^ = 0) Maxwell's equations become 

D'^Ia =2lUj''Ia, (21) 

+ i Cmlla = - |0Xa + (^ab^^ - [uJ, I]a - i [u, T]a ■ (22) 

It follows that the source-free Maxwell equations in an arbitrary curved spacetime, 
relative to an arbitrary congruence of observers, are invariant under the covariant global 
spatial duality rotation Xa c^'^Xq, where is constant. The energy density and 
Poynting vector 

U = = liEaE- + HaH-) , (23) 

Pa=^^^,I]a = [E,H]a, (24) 

are natural group invariants. Their invariance also follows from the duality invariance 
of the energy-momentum tensor p2|, 



Ma'='^{FacF''+ *Fac*F''^) , (25) 

since U = MabU°'u^ and Pa = —M(^a)bu''- Using the identity (5), and the propagation 
equations (19) and (20), we find a covariant energy conservation equation: 

if + D^Pa = - |et/ - 2u''Pa + Oab [E^E'' + H"" H^) . (26) 

This reduces in fiat spacetime for inertial observers to the well-known form dtU+div P = 
0. 

A further natural group invariant is 

TTafe = —I (alb) = —E(aEb) — H(^aHb) , (27) 



which is just the anisotropic electromagnetic pressure It occurs in the last term of 
the conservation equation (26), i.e. — aafevr'**. 
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For later comparison with the gravitational case, we conclude this section 
by considering the propagation of source-free electromagnetic waves on an FLRW 
background, assuming that Ea = = Ha in the background. We linearize and take 
the curl of equation (19), evaluating curl curl i^a by the identity (13) and equation (18). 
We eliminate curlii^a by linearizing equation (20), taking its time derivative, and using 
identity (9). The result is the wave equation 

a^Ha = -Ha + ^^Ha ^ hHHa + [^H^ + \p - p) Ha , (28) 

where p is the background pressure, and we used the FLRW field equation ?>H = 
—3H^ — ^{p + 3p). A similar wave equation may be derived for Ea. 



3. The Bianchi identities and nonlinear duality 

The Maxwell analogy in General Relativity is based on the the correspondence C, 



abed 



Fab, where the Weyl tensor Cabcd is the free gravitational field (see [|T9|)- For a given m", 
it splits irreducibly and covariantly into 

Eab = CacbdU'^u'^ = E(^ab) , Hab = *CacbdU''u'^ = H (^ab) , (29) 

which are called its 'electric' and 'magnetic' parts by analogy with the Maxwell 
decomposition (|T^). These gravito-electric/magnetic spatial tensors are in principle 
physically measurable in the frames of comoving observers, and together they are 
equivalent to the spacetime Weyl tensor, since [^| 



Cab"" = 4 + + 2eabeU^'H'''^- + 2u^aHb]ee"'' . (30) 

This is the gravito-electromagnetic version of the expression (p!6D. The electromagnetic 
interpretation of Eab and Hab is reinforced by the fact that these fields covariantly (and 
gauge-invariantly) describe gravitational waves on an FLRW background (including the 
special case of a fiat vacuum background) |]^, ^ . 

In the 1 + 3 covariant approach to General Relativity , the fundamental quantities 
are not the metric (which in itself does not provide a covariant description), but the 
kinematic quantities of the fluid, its energy density p and pressure p, and the gravito- 
electric/magnetic tensors. The fundamental equations governing these quantities are the 
Bianchi identities and the Ricci identities for m", with Einstein's equations incorporated 
via the algebraic definition of the Ricci tensor Rab in terms of the energy-momentum 
tensor Tab- We assume that the source of the gravitational field is a perfect fluid (the 
generalization to imperfect fluids is straightforward). The Bianchi identities are 

^'^Cabcd = V[a (— -Rfcjc + \Rgb]c^ , (31) 

where R = Ra"^ and Rab = Tab — \Tfgab- The contraction of (|31| ) implies the conservation 
equations. The tracefree part of ( PT| ) gives the gravitational equivalents of the Maxwell 
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equations (17)-(20), via a covariant 1 + 3 decomposition In our notation, these 

take the simphfied form: 

D'^afe = - Suj'Hab + iD,p + [a, H]a , (32) 

D'Hab =3u;'E,b + {p + p)u;a-W,E]a, (33) 

E{ab) — CUll Hab = — QEab + 3(Tc(a-E'b) — Uj" E cd{aEb)^ 

+ 2u'e,d(aHb)'^ - |(P + P)(^ab , (34) 

H{ab) + curl Eab = - QHab + ^CfciaHb)" - U^EcdiaHb)^ 

-2M%,,(,Efe)^. (35) 

These are the fully nonlinear equations in covariant form, and the analogy with the 
Maxwell equations (17)-(20) is made strikingly apparent in our formalism. 

Vorticity couples to the fields to produce source terms in both cases, but gravity 
has additional sources from a tensor coupling of the shear to the field. The analogue 
of the charge density ^ as a source for the electric field, is the energy density spatial 
gradient D^p as a source for the gravito-electric field. Since 'Dap covariantly describes 
inhomogeneity in the fluid, this is consistent with the fact that the gravito-electric field 
is the generalization of the Newtonian tidal tensor 0. 

There is no magnetic charge source for if^, but the gravito- magnetic field Hab has 
the source (p + p)uja. Since p + p is the relativistic inertial mass-energy density 0, 
(p + p)uja is the 'angular momentum density', which we identify as a gravito- magnetic 
'charge' density. Note however that angular momentum density does not always generate 
a gravito-magnetic field. The Godel solution provides a counter-example, where 
Hab = and the non-zero angular momentum density is exactly balanced by the 
vorticity/ gravito-electric coupling in equation (33), with aab = 0. 

For both electromagnetism and gravity, the propagation of the fields is determined 
by the spatial curls, together with a coupling of the expansion, shear, vorticity, and 
acceleration to the fields. The analogue of the electric current ja is the gravito-electric 
'current' (p + p)<Jab, which is the 'density of the rate- of- distortion energy' of the fluid. 
There is no magnetic current in either case. 

If the Maxwell field is source-free, i.e. g = = ja, and the gravitational field is 
source-free, i.e. p = = p, then the similarity of the two sets of equations is even 
more apparent, and only the tensor shear coupling in the case of gravity lacks a direct 
electromagnetic analogue. (Note that these shear coupling terms govern the possibility 
of simultaneous diagonalization of the shear and Eab, Hab in tetrad formulations of 
general relativity |^T], 0.) 



To obtain the gravitational analogue of the complex equations (21) and (22), which 
lead to the Maxwell duality invariance, we consider the vacuum case p = = p. 
In general, u"" is no longer uniquely defined in vacuum, although in particular cases 
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(such as stationary spacetimes), there may be a physically unique choice. However, our 
results hold for an arbitrary covariant choice of m", without any special conditions on the 
congruence. By analogy with the complex electromagnetic spatial vector la, we define 
the complex gravito-electromagnetic spatial tensor 

^ab = Eab + i Hab ■ (36) 

Then equations (32)-(35) reduce to: 

B^Iab = 3i u^Xab - i [o-, T]a , (37) 

Xl^ab) + i curl J„b = - Qlab + ^CTcl^alb)" 

— ^''^cdiaXb)'^ — 2i U^Ecdia'^b)'^ ■ (38) 

Apart from the increased economy, the system is now clearly seen to be invariant under 
the global U{1) transformation: 

lab ^ e^'^Iab, (39) 

which is precisely the tensor (spin-2) version of the vector symmetry of the source- 
free Maxwell equations. We have thus established the existence of the covariant 
spatial duality at the level of the physically relevant gravito-electric/magnetic fields, 
in the general (non-perturbative, arbitrary observer congruence) vacuum case. (As with 
electromagnetism, duality invariance breaks down in the presence of sources.) 

A covariant super-energy density and super-Poynting vector arise naturally as 
invariants under spatial duality rotation, in direct analogy with the Maxwell invariants 
of equations (23) and (24): 

U = \X'''Xab = \ {EabE^' + HabH-') , (40) 
Pa = Y^,Aa = [E, H]a = SabcE^H'^'' . (41) 

This reflects the duality invariance of the Bel-Robinson tensor [|l| 

Mab"" = \ {CaebfC''"^ + *Caebf V''^) , (42) 

which is the natural covariant definition of the super-energy-momentum tensor for the 
free gravitational field, since |10 



U =MabcdU''uW, (43) 
Pa= - M{a)bcduW . (44) 

The agreement between equations (43) and (40) follows obviously from equation (^) 
on using equation (^). However, it is not obvious that equation (44) agrees with our 
equation (41) for the super-Poynting vector, and one requires the identity ( pUj ) to show 
the agreement. 
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Our expression (40) for the gravitational super-energy density gives a direct and 
clear analogy with the electromagnetic energy density (23). Our expression (41) for 
the gravitational super-Poynting vector, in terms of the tensor generalization of the 
vector product, provides a clearer analogy with the electromagnetic Poynting vector (24). 
The analogy is reinforced by the fact that U and Pa obey a super-energy conservation 
equation which is the tensor version of the electromagnetic energy conservation equation 
(26). To show this, we need the new covariant identity (6). Using this and the Bianchi 
propagation equations (34) and (35), we find that 

U + D'^P, = -2eU- Ail-Pa + 3a'^a [^6)c^'^' + H.^.H'^''] . (45) 

This is the non-perturbative and covariant generalization of Bel's linearized conservation 
equation dtU = -divP. 

The last term in the conservation equation (45) contains another natural group 
invariant 

which we interpret as the anisotropic super-pressure of the gravito-electromagnetic field. 
A further covariant quantity that may be naturally constructed from equation (^) 

is 

T TO-b / 771 rpab tt Tjcib\ , <-> ■ 771 Tjab 

-Lab-L — [J^ab-l^ —J^abJ^ ) + Zl H^abJ^ 

= liCabcdC^"''' + iCabcdC-'"") (47) 

which is not invariant under equation (^). However it proves very useful in categorizing 
spacetimes and vanishes in Petrov type-Ill and type-N spacetimes [T^ (supporting 



the existence of gravitational waves in these spacetimes, since the analogous quantities 
vanish for purely radiative electromagnetic fields). Further, the electromagnetic 
analogue of the real part of equation (47), namely EaE'^ — HaH"^ = —\FabF°'^ , is 
just the Lagrangian density. The analogue of the imaginary part is EaH"^ = jFab*F°-^ 
whose integral in non-Abelian gauge theories is proportional to the topological instanton 
number. 

As pointed out in the introduction, duality rotations preserve the Bianchi identities 
in vacuum, but not the Ricci identities for u"". This is clearly apparent from the spatial 
tensor parts of the Ricci identities 0, which in our formalism have the simplified form 

Eab = D(a^b) - 0-(a6) - ^Qc^ab " (^ciaC^b)"" " (^{a^Jb) + U{aUb) , (48) 

Hab = curl aab + D^^u;;,) + 2u(^a^b) ■ (49) 

In order to preserve the Ricci identities, and map Einstein solutions to Einstein solutions, 
one needs to perform kinematic transformations in addition to the duality rotation. An 
example is presented in the following section. 
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The electromagnetic analogy suggests a further interesting interpretation of the 
kinematic quantities arising from the Ricci equations (48) and (49). f In flat spacetime, 
relative to inertial observers, the electric and magnetic vectors may be written as 

E = VV — dtd , H = curl a , 

where V is the electric scalar potential and a is the magnetic vector potential. f 
Comparing now with the Ricci equations (48) and (49), we see that the four-acceleration 
is a covariant gravito-electric vector potential and the shear is a covariant gravito- 
magnetic tensor potential. The vorticity derivative in (49) has no electromagnetic 
analogue, and vorticity appears to be an additional gravito- magnetic vector potential. 
Furthermore, the gauge freedom in the electromagnetic potentials does not have a direct 
gravitational analogue in the Ricci gravito-potential equations (48) and (49), since the 
gravito-electric/magnetic potentials are invariantly defined kinematic quantities. (Note 
that the Lanczos potential for the Weyl tensor does have a gauge freedom analogous to 
that in the Maxwell four-potential |]12[.) 



The remaining Ricci equations in 1 + 3 covariant form are |19 



+ 1©' = - |(P + 3p) + B^Ua + u^Ua + 2a;'^cu, - a'^V^b , (50) 

UJl^a) + I'duJa = - |curl Ua + (TabUJ^ , (51) 

fD^e = -cm\uJa + B^aab + 2[uj,u]a, (52) 

D'^C^, = U^'Ua , (53) 

and do not involve the gravito-electromagnetic field. 

Finally in this section, we extend the analogy to wave propagation. The magnetic 
wave equation (^) has a simple gravito-magnetic analogue. In order to isolate the purely 
tensor perturbations of an FLRW background in a covariant (and gauge-invariant) way, 
one imposes Ua = |T^. We hnearize and take the curl of equation (34), using the 
linearizations of equations (35) and (33), and identities (10) and (14). This does not 
directly produce a wave equation, since the curl of the shear term in (34) has to be 
ehminated. (In the Maxwell case this feature did not arise, since we set ja = 0.) The 
elimination is achieved via the Ricci equation (49), and we find that 

n^H^, = -Hab + B'H^t, ^ 7HH,t, + 2 (SH' - p) H^b , (54) 



in agreement with |[T^, |T8[ , and in striking analogy with the magnetic wave equation (|28 



Further discussion of covariant gravitational wave theory may be found in |T^, |T8|, |T9|, |43[] 



f Note that these Ricci equations have the same form in the non-vacuum case. 

f The covariant form of these potentials is = u'^Aa, eta — ^(a>j where Aa is the four-potential. 
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4. Gravitational monopoles 

The electromagnetic correspondence we have developed suggests a covariant 
characterization of gravito-electric (magnetic) monopoles, as stationary vacuum 
spacetimes outside isolated sources, with purely electric (magnetic) free gravitational 
field, i.e., Hab = {Eab = 0). This is reinforced by the fact that monopoles do not 
radiate, and gravitational radiation necessarily involves both Eab and Hab nonzero (see 
T8| , p!9| , consistent with Bel's criterion 7^ [I], |10|)- Our identification in the 



previous section of density inhomogeneity and angular momentum density as sources 
of, respectively, gravito-electric and gravito-magnetic fields, suggests that the monopole 
sources will be respectively mass and angular momentum. However, as pointed out 
previously, it is possible that non-zero angular momentum is compatible with a purely 
gravito-electric field, as illustrated by the Godel solution. 

The four-velocity field is not defined by a fluid, but is defined as the normalization 
of the stationary Killing vector field = As a consequence of Killing's equations. 



we have = = aab [0, so that 



VftMa = Eabc^^" - UaUb ■ 

The covariant equations governing non-perturbative monopoles are complicated. Some 
simplification arises from the Killing symmetry, which implies 



Ci:UJa = iiOa + UaUj'Bbi = , 

^^Hab = iHab + 2^Uj''ecd{aHhf — 2^U(^aHb) 

and a similar equation for Eab- Then it follows that 



LOl^a) = , (55) 

H{ab) = — 2uj''ecd(aHb)'^ , (56) 

E{ab) = — 2uj''ecd{aEbf ■ (57) 

Now equations (55)-(57), together with the basic monopole conditions, are apphed 

to the Bianchi equations (32)-(35) and Ricci equations (48)-(53). We obtain: 

Gravito-electric and -magnetic monopoles: 

B'^ila = -ifUa-2uo''uOa. (58) 

Cj{a) = , (59) 

curl Ua = , (60) 

CUlluJa = - 2[u,Uj]a, (61) 

D'^Ua =u''iOa. (62) 
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Gravito- electric monopole: 

B'E,, =0, (63) 

=Ea,u;\ (64) 

E^ab) = , (65) 

= u'EcdiaEb)'' , (66) 

curlE^fc = -2n%,,(,Efe)^ (67) 

Eab - D(aMb> = U{aUb) " ^^{a^b) , (68) 

^{a^b) = - "iu^ai^b) ■ (69) 

Gmvito-magnetic monopole: 

D'Hab =0, (70) 

=H,bio\ (71) 

H(^ab) =0, (72) 

= iu'^e,d^aHb)\ (73) 

curli7„6 = -2u'ecdiaHb)\ (74) 

D(aMb) = - M(aMfe> + i^{a^b) , (75) 

-f^afe - D(at^b) = '^U{a^b) ■ (76) 

Equation (60) implies that there exists an acceleration potential: 

Ua = . (77) 



This holds even when Ua ^ 0, despite the identity (4), since $ is invariant under 
so that $ = 0. Equation (61) shows that curlci;a is orthogonal to the vorticity and 
four-acceleration: 

CiJ^CUrl UJa = = M^curl Ua . 

Schwarzschild spacetime, where also Ua = (since staticity implies u"" is hyper- 
surface orthogonal), is clearly a non-perturbative gravito-electric monopole according 
to our covariant definition: it is a static vacuum spacetime satisfying Hab = 0, by virtue 
of the Ricci equation (49). Equations (|7^) and (58) imply 

+ D''<1'D„<1) = . (78) 

The solution $ determines Ua and Eab, and equation ( [7^ ) ensures that the monopole 
conditions (58)-(69) are identically satisfied. 

It is not clear whether there exist consistent non-perturbative gravito-magnetic 
monopoles, i.e. spacetimes satisfying the covariant equations (58)-(62) and (70)- 
(76). t However, linearized gravito-magnetic monopoles have been found, for example 

t In it is shown that non-flat vacuum solutions with purely magnetic Weyl tensor are a very 
restricted class, and it is suggested that there may be no such solutions. 
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the Demianski- Newman solution (see below). It is also not clear whether there exist 
gravito-electric monopoles with angular momentum (i.e. Ua ^ 0). 

In the case of linearization about a flat Minkowski spacetime, the right-hand sides 
of equations (58)-(76) may all be set to zero. In particular, equation (61) implies that 
there is a vorticity potential: 

COa ~ D,^ . (79) 

The linearization of equations (68) and (76), together with the scalar potential equations 
{\n\) and ([79|), then imply that the curls vanish to linear order. Thus the linearized 
gravito-electric monopole is covariantly characterized by equations (|77D, (ffOD and 

D^^^O, Eat-B^D,^, D^,D,)^^0, (80) 

while for the linearized gravito-magnetic monopole 

D^^^O, Hab^DaD,^, D^,Db)$^0. (81) 

It follows in particular that a linearized non-rotating gravito-electric monopole is 
mapped to a linearized non-accelerating gravito-magnetic monopole via 

lab -^ilab, i^a-^Ua, Ua ^ "^^a • (82) 

Linearized Schwarzschild spacetime is readily seen to satisfy equation (80) with $ = 
—M/r, where M is the mass and r the area coordinate. Using the spatial duality rotation 



and kinematic interchange described by equation (|82D , this monopole is mapped to a 
linearized non-accelerating gravito-magnetic monopole with potential = —M/r. In 
comoving stationary coordinates, the metric of the linearized magnetic monopole follows 
from iia = and Ua = Da\l/, using a theorem in (p 24): 

ds^ = -df + dr^ + (de'^ + sin^ 6 dip^) + 4M cos Odipdt. (83) 

This is a Taub-NUT solution with m = 0, £ = —M and linearized in I ([0], p. 133; 
see also ||5B|). In fact, this is precisely the linearized solution found in so that 



we have a covariant characterization of that solution in the framework of gravitational 
duality. Clearly the magnetic 'charge' M is an angular momentum parameter, not 
a mass parameter, and the metric in equation (83) describes an isolated source with 
angular momentum but no mass. 

5. Concluding remarks 



A covariant 1 + 3 approach, based on ^ and its extension |T9|], is ideally suited 
to an analysis of the free gravitational field that is based on observable physical 
and geometric quantities, with a clear and transparent analogy in well-established 
electromagnetic theory. We have used such an approach, including in particular the 
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generalization of covariant spatial vector analysis to spatial tensor analysis, which 
involves developing a consistent covariant definition of the tensor curl and its properties. 
Via this approach, we showed the remarkably close analogy between the Maxwell 
equations for the electric/magnetic fields and the Bianchi identities for the gravito- 
electric/magnetic fields. Although this analogy has long been known in general terms, 
our approach reveals its properties at a physically transparent level, with a detailed 
accounting for each physical and geometric quantity. We found new interpretations of 
the role of the kinematic quantities - expansion, acceleration, vorticity and shear - in the 
source and coupling terms of gravito-electromagnetism. The tracefree part of the Ricci 
identities also reveals the role of the kinematic quantities as gravito-electric/magnetic 
potentials. 

The analogy provides a simple interpretation of the super-energy density and 
super-Poynting vector as natural U{1) invariants, and we derived the exact nonlinear 
conservation equation that governs these quantities, and which involves a further natural 
invariant, i.e. the anisotropic super-pressure. We also used the analogy to show 
that a covariant spatial duality invariance exists in vacuum gravito-electromagnetism, 
precisely as in source-free electromagnetism. Duality invariance has been important in 
some recent developments in field and string theory, and the gravito-electromagnetic 
invariance in the form found here may also facilitate new insights into gravity. A 
crucial feature in the gravitational case, arising from its intrinsic nonlinearity, is that 
the duality invariance does not map Einstein solutions to Einstein sohitions, since the 
Ricci identities are not invariant. Further work is needed to investigate whether a 
simultaneous geometric or kinematic transformation can be found, so that the Bianchi 
and Ricci equations are invariant under the combined transformation. 

We showed that in linearized vacuum gravity, there is a simple combined duality/ 
kinematic transformation that maps the Schwarzschild gravito-electric monopole to 
the Demianski-Newman gravito-magnetic monopole. This covariant characterization 
of the relation between these linearized solutions was based on our covariant definition 
of gravito-monopoles in the general nonlinear theory. Further work is needed on the 
governing equations for these monopoles, in particular to see whether nonlinear gravito- 
magnetic monopole solutions may be found. A better understanding of the relation 
between nonlinear gravito-electric/magnetic monopoles could, as in field theory, open 
up new approaches and insights. 
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